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ABSTRACT 


The Building and Fire Research Laboratory (BFRL) was given the opportunity to make measurements 
during fire calibration tests of the heat detection system in an aircraft hangar with a nominal 30.4 m 
(100 ft) ceiling height near Dallas, TX. Fire gas temperatures resulting from an approximately 
8250 kW isopropyl alcohol pool fire were measured above the fire and along the ceiling. The results 
of the experiments were then compared to predictions from the computer fire models DETACT-QS, 
FPETOOL and LAVENT. In section A of the analysis conducted, DETACT-QS and FPETOOL 
significantly underpredicted the gas temperatures. LAVENT at the position below the ceiling 
corresponding to maximum temperature and velocity provided better agreement with the data. For 
large spaces, hot gas transport time and an improved fire plume dynamics model should be 
incorporated into the computer fire model activation routines. A computational fluid dynamics (CFD) 
model, HARWELL FLOW3D, was then used to model the hot gas movement in the space. 

Reasonable agreement was found between the temperatures predicted from the CFD calculations and 
the temperatures measured in the aircraft hangar. In section B, an existing NASA high bay space was 
modelled using the computational fluid dynamics model. The NASA space was a clean room, 27.4 m 
(90 ft) high with forced horizontal laminar flow. The purpose of this analysis is to determine how the 
existing fire detection devices would respond to various size fires in the space. The analysis was 
conducted for 32 MW, 400 kW, and 40 kW fires. 


Key words: clean rooms; computer models; computational fluid dynamics; detector response; field 
modeling; fire detection; fire models; forced air flow; fire plumes; fire tests; high bays; response 
time; sprinkler response 
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1.0 INTRODUCTION 


Large spaces, such as those found in warehouses, historical buildings, atriums, and aircraft hangars, 
represent some of the most difficult fire protection challenges since they are frequently of historical 
significance, contain large quantities of fiiel, have unique geometries, and/or present special life safety 
problems. Accurate detector activation predictions are important in these large spaces, as timely 
detection of a fire is more difficult due to the distance heat and products of combustion must travel to 
reach sprinklers and detectors.^. Since fires frequently grow at an exponential rate, an increased time 
to detection results in larger fires at the time of detection and larger fires to be suppressed by, e.g., 
an automatic sprinkler system. Even a modest uncertainty in the prediction of the activation time may 
lead to a large uncertainty in the fire size used to predict the hazard, a central element of many fire 
protection analyses. 

There has been substantial effort to verify activation predictions in small and medium sized rooms, 
but little in large spaces. In January, 1993, a report describing model verification tests conducted in a 
14 m (46 ft) high hangar was released by NIST [1]. Although the correlations in the computer fire 
models were based on ceiling heights up to 16 m (52 ft), previous to this study, comparisons of 
computer fire model predictions with experimental measurements for ceiling heights approaching that 
height had not been made. Conducting verification experiments for large spaces is difficult due to the 
lack of availability of adequate facilities for live fire tests. The largest dedicated fire test facility in 
the United States is limited to a height of 18.2 m (60 ft). 

Why is model verification important? It is the responsibility of the user of a model to determine the 
suitability of a particular fire model for a given situation. A successful comparison of the model 
predictions with real-scale fire results helps the user make that decision. Comparison of a fire model 
against appropriately chosen experiments can provide the necessary base to extend the use of the 
model to real fire applications. Comparisons also allow model developers to examine the accuracy of 
a model and identify the needs for future model development. It is critical that the user have 
comparisons of model predictions to a wide range of experimental data. 

The Building and Fire Research Laboratory (BFRL) was given the opportunity to make measurements 
during fire calibration tests of the heat detection system in an aircraft hangar with a nominal 30.4 m 
(100 ft) ceiling height near Dallas, TX. Three closed-door tests were conducted 

This project was funded in two phases. The U.S. General Services Administration (GSA) provided 
initial funding to conduct the experiments and to compare measured data with the predictions of 
DETACT-QS, LAVENT, and FPETOOL computer fire models. GSA’s interest is to verify and 
provide recommendations for the use of the GSA Engineering Fire Assessment System for predicting 
the response of sprinklers and detectors in large spaces. FPETOOL is the main component of this 
system. The National Aeronautics and Space Administration (NASA) provided funding to compare 
the predictions of a computational fluid dynamics, or field model to the experiment. Comparisons of 
the experimental data to both the computer fire models and the field model are described in section A 
of this report. NASA also funded section B of this project which uses the field model to predict 
smoke movement in a NASA Goddard high bay space. NASA’s interest is in the ability of field 
modelling to provide information that would be useful in selecting effective fire protection options for 
their high bay spaces. 
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2.0 EXPERIMENTAL 


2. 1 Building layout 

The aircraft hangar, pictured in figure 1, measured approximately 389 m (1280 ft) long by 115 m 
(377 ft) deep (81 m (270 ft) to a firewall), and 30.4 m (100 ft) high. The hangar had seven bays and 
was capable of housing seven wide-bodied aircraft side by side. The hangar contained draft curtains 
which were spaced approximately 12.5 m (41 ft) apart as shown in figure 2. Draft curtains are 
partitions extending vertically down from and perpendicular to the ceiling. Draft curtains are 
designed to slow the spread of combustion gases across the ceiling. The hot gases collect between the 
draft curtains and flow downward until spilling into the adjacent curtained areas. This concentration 
of smoke and hot gases inside the draft curtain is intended to speed up the operation of sprinklers 
within the curtained area. Draft curtains in this facility were 3.7 m (12 ft) deep. A plan view of the 
fire test curtained area and the two adjacent curtained areas is shown in figure 2. 

The ceiling geometry is shown in figure 3. Perpendicular to the draft curtains, a series of large I- 
beams are located near the ceiling which, when coupled with the draft curtain geometry, appear to 
partition the curtained areas into cubes. A second set of I-beams runs both below the draft curtains 
and the I-beams forming a second set of open cubes at a height of 24 m (79 ft) off the floor. Large I- 
beams join the top and bottom I-beams in a triangular structure. The underside of the roof is made of 
corrugated steel with 38 mm (1.5 in) deep corrugations. 

The front of the hangar consisted primarily of sliding doors, which if all were open would provide a 
maximum opening 289 m (948 ft) in length and 22 m (72 ft) high. The hangar also contained many 
personnel doors both to the outside and to other parts of the building. The large doors were closed 
during the tests; however, some of these personnel doors may well have been open during the tests. 
(There was much work going on in the building during testing.) The building HVAC system was off 
for ftie duration of the three tests. 

2.2 Test fire 

The test fire was an array of nine pans, each 0.91 m x 0.91 m (3 ft x 3 ft), placed together to form a 
square fire area totalling 7.5 m 2 (81 ft 2 ) as shown in figure 4. The effective circular diameter of this 
fire is 3.1 m (10 ft). The pans were 0.076 m (0.25 ft) deep. The pans were supported on bricks, 

127 mm (0.42 ft) off the floor. The test fuel was technical grade isopropyl alcohol. A total of 135 
liters (35.1 gal) of fuel was used per test, 15 liters (4.0 gal) of fuel per pan. The fuel depth in each 
pan was approximately 13 mm (0.5 in). Cloth wicks were draped between the pans to aid in fire 
spread after ignition of one pan. 

The fire was located on the floor in the center of the hangar building. At 8 seconds after ignition of 
the first pan, all the fuel pans were fully involved. A picture of the fire at steady-state is shown in 
figure 5. Steady burning was maintained for approximately three and one-half minutes, after which 
the flame height started to decrease. The fire was allowed to bum out, which took several additional 
minutes. 

As there was no load cell or oxygen calorimeter available for use dining these experiments, the heat 
release rate of the fire was estimated from previous bum tests conducted at NIST. These tests are 
described in [1]. Isopropyl alcohol pan fires 3.3 m 2 (36 ft 2 ) were conducted under the NIST large 
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calorimeter. These tests used 4-0.91 x 0.91 m (3 x 3 ft) steel pans. The effective circular diameter 
of this fire was 2.1 m (6.9 ft). In the bum tests conducted, die average burning rate was 0.036 
kg/m 2 /s (0.0074 lb/ftVs), and the average heat release rate was 3650 kW (3460 Btu/s). The net heat 
of combustion is reported in the NFPA handbook [2] as 30.45 MJ/kg (13107 Btu/lb). 

Assuming there is sufficient oxygen available for complete combustion, the heat release rate of a fire 
can be calculated as: 

0 = rh" * Ah c * A f 

where: 

0 — heat release rate (kW) 
m" = burning rate (kg/m 2 /s) 

Ah c = net heat of combustion (kJ/kg) 

A f = surface area of fuel (m 2 ) 

To characterize the burning rate of a liquid pool fire, it is required that a classification be made 
according to the dominant heat transfer mechanism, which can be dependent on pool diameter [3]. 
When the pool diameter, D, is very small, D < 0.03 m (0.01 ft), conductive heat transfer determines 
the rate of burning, while the radiative term predominates if D is large, D ^ 1.0 m (3.3 ft). At D = 
1.0 m (3.3 ft) and above, the burning rate has been shown to be constant [4]. Since the isopropyl 
alcohol pan fires conducted at NIST were for an effective diameter of 2.1 m (6.9 ft), which is well 
within the radiative dominated range, the same burning rate was used for the 3.1 m (10 ft) effective 
circular diameter fire conducted in the aircraft hangar. Using the above equation, the calculated heat 
release rate for the fire in the aircraft hangar is approximately 8250 kW (7820 Btu/s). 

The heat release rate as a function of time is shown in table 1. An initial growth phase of the fire 
was selected to provide a reasonable fit with the measured temperature rise early in the fire, as well 
as with test observations. 

The fraction of the heat release rate that is radiated from the flaming region directly affects the gas 
temperature. The radiant heat fraction is sometimes referred to as X,.. For many common fuels with 
lu min ous flames, it has been found that the radiative fraction is 35 percent [5]. Although many 
alcohols bum without luminous flames, isopropyl alcohol flames are very luminous, appearing much 
like the flames produced from burning heptane. Since isopropyl alcohol bums with luminous flames, 
the analysis was performed with a radiative fraction of 35 percent. 

To determine the sensitivity of the gas temperature predictions to heat release rate, a 20 percent 
higher heat release rate, 9900 kW (9375 Btu/s) was also examined in this report. This provides a 
higher end prediction of the temperatures in the hangar. The higher heat release rate analysis takes 
into account possible uncertainties in measurement of prior tests of the burning rate performed in the 
NIST large calorimeter. 
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2.3 Measurements 


2.3.1 Type of measurements taken 

Thermocouples were used to measure the gas temperatures above the fire and along the ceiling. The 
thermocouples were b are-bead, 0.52 mm (0.02 in) diameter, type K, chromel-alumel thermocouples. 
Brass disks were used to simulate the heating of fusible links. The disks were 9.5 mm (0.38 in) 
diameter, 0.81 mm (0.03 in) thick brass. The RTI of the brass disks is 33.4(m-s) 1/2 (60.5 (ft-s) ,/2 ). 
Figure 6 shows the instrumentation layout. The temperatures of the thermocouples and brass disks 
were recorded every 4 seconds on a computer-based data acquisition system. 

To measure the gas temperature , six temperature measurements ware taken vertically below the 
ceiling east of the fire at radial distances of 1.5, 3.0, 4.6, and 6.1 m (5, 10, 15, and 20 ft) from the 
centerline of the fire and west of the fire at a radial distance of 3.0 m from the centerline of the fire. 
At each radial position, a thermocouple and calibrated brass disk were placed at 0. 15 m (0.5 ft) and 
0.91 m (3 ft) below the ceiling. Thermocouples were also placed at 0.30 m (1.0 ft) and 1.5 m (5 ft) 
below the ceiling. Assuming that the smoke and temperature flow together, the thermocouple 
readings can indicate the initial arrival of the smoke to the thermocouple location. This assumption is 
supported by the experiments of Heskestad and Delichatsios [6J. 

To track smoke filling in the fire test curtained area, thermocouples wore placed to the east and west 
of the fire at 5.5 m (18 ft) radially out from the fire centerline. Thermocouples were placed 
vertically from the ceiling to the bottom of the draft curtain as shown in figure 6 to determine the 
layer height. 

To measure the centerline plume temperatures, eight thermocouples were placed directly over the fire 
at 0.025, 3.0, 6.1, 9.1, 12.2, 15.2, 18.2, and 21.3 m (0.082, 10, 20, 30, 40, 50, 60, and 70 ft) 
below the ceiling. 

To track the flow of smoke, a thermocouple was placed 0.025 m (0.082 ft) directly under the draft 
curtains between the test area and the adjacent curtained areas. An array of four thermocouples was 
placed at the center of the east adjacent curtained area and of the west adjacent curtained area. 
Thermocouples were also placed 0.025 m (0.082 ft) below the draft curtains at the far end of the east 
adjacent curtained area and at the far end of the west adjacent curtained area. Thermocouples were 
placed to the north and south of the fire, at radial distances of 3.0, 6.1, and 9.1 m from the pan 
centerline (10, 20, and 30 ft), each 0. 15 m (0.5 ft) below the ceilin g 

2.3.2 Uncertainty of measurements 

The manufacturer of the thermocouple wire [7] states that there is a ± 2.2 °C uncertainty in the 
thermocouple reading. The resolution of the data acquisition system is 0.3 °C. Experimental 
uncertainties such as re-radiation and chemical dissociation are negligible for the temperature range 
realized in this experiment. 

The fire pulsations are reflected in the temperature fluctuations seen in the curves for measured 
temperature. A mean temperature for each measurement point was calculated and is reported with the 
associated standard deviation. This deviation is plus or minus three to seven degrees Celsius 
depending on the location of the thermocouple. 
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3.0 COMPUTER MODELS 


The NFPA handbook groups computer fire models into two categories: enclosure fire models and 
special-purpose models [8]. Enclosure fire models are structured to predict the interaction of multiple 
fire processes involving heat transfer, fluid mechanics, and combustion chemistry occurring in an 
enclosure. Special purpose computer models are computer models designed for special purpose 
analyses such as: structural fire resistance, prediction of response time of heat detectors and automatic 
sprinklers, egress design, or performance of smoke control or ventilation systems. Fire models are 
farther subdivided into probabilistic and deterministic models. Probabilistic models treat fire growth 
as a series of sequential events or states. Probabilities are assigned to each transfer point based on 
analysis of relevant experimental data and historical fire incident data. In contrast, deterministic 
models represent the processes encountered in a compartment fire by interrelated mathematical 
expressions based on physics and chemistry. According to this system classification, DETACT-QS, 
LAVENT, and the activation subroutine in FPETOOL are all special purpose deterministic fire 
models [8]. 

A CFD model can be used as a deterministic enclosure model. A CFD model solves the mass, 
energy, and momentum equations directly which provides more flexibility in modelling complex 
situations. In using a CFD model, the enclosure is divided up into many small spaces or grid 
volumes and the conservation equations are then solved to yield the values of temperature, velocity, 
etc. at the center of these grid volumes. CFD modelling then provides temperature and velocity 
distributions across the enclosure. 

3.1 DETACT-QS and FPETOOL 

It is important to understand the physics of the computer models and the assumptions built into each 
code. FPETOOL, written by Nelson [9,10] and DETACT-QS (DETector ACTuation Quasi-Steady), 
written by Evans and Stroup [11,12], are each based on experimental correlations developed by 
Alpert for steady-state fires [13]. These correlations give die maximum temperature and maximum 
velocity of the ceiling jet as a function of the steady state heat release rate of the fire, the radial 
distance from the fire to the detector, and the height of the ceiling above the fuel. Both programs 
assume the detector or sprinkler is located at a distance below die ceiling corresponding to maximum 
temperature and velocity. Alpert’s correlations assume a smooth, unconfined ceiling. They also 
assume that steady-state correlations can be applied to a growing fire over small time intervals. Both 
DETACT-QS and FPETOOL perform a linear interpolation to determine the heat release rate at 
intermediate times. 

In both programs, the transport time of the smoke and hot gases from the fire to the thermal detector 
is neglected. Also in both programs, the detector is subject to the maximum temperature and velocity 
of the ceiling jet. FPETOOL accounts for the impact of the hot gases entrained into the ceiling jet on 
the temperature and velocity of the jet as it passes through the hot smoke layer; DETACT-QS does 
not. The expressions used in DETACT-QS and FPETOOL for temperature and velocity are based on 
a nondimensional parameter r/h, where r is the radial distance from die centerline of the fire to the 
thermal device, and h is the height from the base of the fire to the ceiling. These expressions are 
independent of radius for r/h £ 0.18, which is assumed to be within the plume region. 

The user inputs for DETACT-QS are shown in table 2. The user inputs for FPETOOL are shown in 
table 3. FPETOOL version 3.01 was used. In both programs, 999 °C is used for the response 
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temperature of the device as the predicted thermal device temperature was compared to the measured 
temperature of the disk and not to the activation time of a sprinkler or detector, i.e., the device was 
programmed not to go off. FPETOOL and DETACT-QS do not model draft curtains. DETACT-QS 
has no walls and assumes open burning. It can not accommodate heat losses to the ceiling. To model 
a draft curtain using FPETOOL, the dimensions of the curtained area are used for the room spacial 
dimensions, the program assumes there are walls. The door/window option is then used to simulate 
the area under the draft curtains as large openings in the walls. 

3.2 LAVENT 

LA VENT (Link-Actuated VENT), is a program written by Davis and Cooper [14,15], to simulate the 
environment and response of sprinkler links in compartment fires with draft curtains and fus'ble-1 ink- 
actuated ceiling vents. LAVENT is similar to FPETOOL and DETACT-QS in that it assume that 
steady state correlations can be applied to a growing fire over small time intervals; and it too neglects 
the transport time of the smoke and hot gases from fire to thermal detector. LAVENT does account 
for the impact of the hot upper layer on the ceiling jet. The important difference between LAVENT 
and the other activation models, is that LAVENT accounts for position of the thermal detector below 
the ceiling in the ceiling jet. It cannot model recessed or concealed sprinklers, as the velocity in 
LAVENT is zero at the ceiling. The expressions for ceiling jet temperature and velocity used in 
LAVENT are independent of radius for r/h £ 0.2, which is assumed to be within the plume region. 
The expression for ceiling jet temperature and velocity as a function of height beneath the ceiling used 
in LAVENT are not valid in the plume region, where there are no appreciable temperature or velocity 
gradients below the ceiling. LAVENT will only provide an accurate temperature prediction in the r/h 
0.2 region at the point below the ceiling where the maximum temperature and velocity occur. This 
unfortunately must be determined by the user through trial and error, although some guidance is 
provided in the user guide. 

LAVENT utilizes draft curtains instead of walls. LAVENT models just one curtained area. When 
that area fills with smoke it spills out into an infinite open space. The user inputs for LAVENT are 
shown in table 4. 

3.3 FLOW3D 

A computational fluid dynamic (CFD) or field model, HARWELL FLOW3D [16], was used to 
simulate the fires in the large spaces analyzed in this report. This model uses a finite difference 
method to solve the Navier Stokes equations in either two or three dimensions on a user-specified 
grid. The Navier-Stokes equations include mass, momentum, and energy conservation. 

While the CFD model enjoys a degree of flexibility that is missing in the special purpose activation 
models, there are a few drawbacks to using a CFD model. The model is typically more difficult to 
use in that an appropriate grid must be chosen to represent the physical space, accurately solve the 
partial differential equations, and provide for adequate representation at the boundaries. These 
problems typically require thousands of grid volumes to provide accurate results. Long computation 
times are therefore required for CFD calculations. Fire-based CFD calculations require a turbulence 
model in order to simulate the gas momentum and energy flow in turbulent regions. The turbulence 
models now available either do not do an adequate job modelling fire plumes or have not been 
thoroughly tested over a large range of fires. 
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For the calculations reported here, the tire is represented by a heat release rate with the energy input 
being distributed over several grid volumes to avoid an unrealistically large energy input in any single 
grid volume. The ideal gas law is used to represent the equation of state for air. Buoyancy is 
incorporated in the calculation by assuming that the gas is compressible and that the density is given 
by the ideal gas law. Fire induced turbulence is modeled using the k-c model 16]. The turbulence 
model depends on a number of user specified parameters which are given in table 5. Wall boundary 
conditions are treated using the standard model for turbulent flow. Table 6 lists the key words used 
in the command tile to activate physical models, the grid volume dimensions for some of the 
calculations, the ambient variables, and the key words to specify the boundary conditions used for the 
FLOW3D calculations. 


4.0 COMPARISON OF COMPUTER FIRE MODELS TO MEASUREMENTS 

4. 1 Plume symmetry 

Table 7 shows the measured mean temperature and standard deviation of the gas at 3.0 m (10 ft) to 
the east of the pan centerline and 3.0 m (10 ft) to the west of the pan centerline at a depth of 0.15 m, 
0.30 m, 0.91 m, and 1.5 m (0.5 ft, 1 ft, 3 ft, and 5 ft) below the ceiling for the three experiments 
conducted. By comparing the east-west temperature measurements, it can be seen that the smoke 
plume was most symmetrical in test 2. Based on plume symmetry, test 2 was used for comparison to 
the predictions of the computer models. 

4.2 Gas temperatures 

Table 8 shows the steady state gas temperatures at 0, 1.5, 3.0, 4.6, 6.1, and 12.5 m (0, 5, 10, 15, 

20, and 41 ft) radial distance from the centerline of the fire as measured in the experimental section. 
The measured temperature reported is the mean temperature from 32 to 100 s. The number in 
parentheses following the mean temperature is the standard deviation. This time period represents 
quasi-steady conditions of the temperature near the ceiling. 

4.2.1 Inside the fire test curtained area 

Figures 7 and 8 show the graphical comparison of the measured and predicted gas temperatures at a 
radial distance of 1.5 m (5 ft) and 6.1 m (20 ft) from the centerline of the fire respectively. These 
radial distances correspond to non-dimensional radii (r/h) of 0.05, and 0.2. The measured gas 
temperatures are shown as solid lines. In each graph, these four solid lines represent the measured 
temperatures at 0.15 m, 0.30 m, 0.91 m, and 1.5 m, (0.5 ft, 1 ft, 3 ft, and 5 ft) below the ceiling. 
Since the temperatures at these four distances below the ceiling are coincident, they were not 
individually labeled. These curves fluctuate, which reflects the pulsation of the fire. The predictions 
of DETACT-QS and FPETOOL are also shown as solid lines. The predicted values are smooth and 
straight. Predictions of gas temperature by LAVENT are shown as various style dashed lines. 

Arrows on the graph point to predictions of LAVENT at vertical distances below the ceiling of d = 
0.15 m, 0.30 m, 0.91 m, and 1.5 m (0.5 ft, 1 ft, 3 ft, and 5 ft). There is only one curve shown for 
DETACT-QS and FPETOOL at each radial distance as these programs do not account for the vertical 
distance below the ceiling. 
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Measurement locations from the fire centerline through 6. 1 m (20 ft) radial distance from the fire are 
within the fire test curtained area as shown in figure 6. For these locations, the r/h values are all 0.2 
or less. For r/h £ 0.2, the temperature predictions of the computer fire models do not change. 
DETACT-QS and FPETOOL underpredict the excess temperature (temperature rise above the ambient 
value of 22 °C) by 34 - 41 percent (depending on the radial distance from the fire centerline) when 
compared to the measured values averaged over a 32 - 100 second time interval. It can also be 
observed from these graphs, that the predictions of LAVENT nearer to the ceiling are closer to the 
measured values, and thus more accurate than the DETACT-QS or FPETOOL values. However, for 
positions further from the ceiling, the predictions of LAVENT are less accurate than DETACT or 
FPETOOL. This is due to the measurement locations being inside the plume. The correlations for 
ceiling jet velocity and temperature as a function of distance below the ceiling used in LAVENT are 
not valid inside the plume, and thus the user must be careful not to use this parameter for locations 
within the plume. This phenomena is illustrated in figure 9. 

However, for high bay areas, often the location of interest is an r/h ^ 0.2 due to the ceiling height. 
Therefore, when running LAVENT inside the plume, the user should use only the prediction at the 
position below the ceiling corresponding to the maximum temperature and velocity. This position 
must be found through trial and error, although some guidance is provided in [14]. LAVENT at the 
point below the ceiling where the temperature is maximum (approximately 0.15 m below the ceiling 
in this case) underpredicts the excess temperature by 6 - 16 percent when compared to the measured 
values averaged over a 32 - 100 second time interval. Future versions of LAVENT should be 
upgraded so that the computer program does not calculate for a ceiling jet (temperature gradient in the 
vertical direction) for r/h ^ 0.2. 

Considering the above mentioned problems with the predictions of the computer fire models, it may 
be useful to use the plume correlations for the plume centerline temperatures to determine the 
temperature rise above ambient for locations inside the plume region. These plume correlations 
assume burning in the open (i.e., no ceiling), and thus free entrainment. There are a number of 
prediction methods available for plume temperature, in general producing very similar results. For 
comparisons performed here, the prediction methods of McCaffrey were selected [18] 



where: 

AT 0 = centerline temperature rise above ambient (K) 

T 0 = ambient temperature (K) 
g = acceleration due to gravity (m/s 2 ) 

C = constant, 0.9 
z = height above base of fire (m) 

Q = heat release rate of die fire (kW) 

k, and ij are constants whose value depends on weather the location of interest is within the flame, 
intermittent, or plume region of the fire. This is determined by die factor z/Q 275 (m/kW 275 ). If z/Q 275 
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is < 0.08, it is considered that this location is in the flaming region of the fire. In this region, k = 
6.8 m 1/2 /s, and rj = 1/2. If z/Q 275 is between 0.08 - 0.2, it is considered that this location is in the 
intermittent region of the fire. In the intermittent region, k = 1.9 m/kW ,/5 s, and ij = 0. If z/Q 275 is 
> 0.2, it is considered that this location is in the plume region of the fire as is the case for the 
calculations near the ceiling. In the plume region of the fire, k * 1.1 m 4/3 /kW 1/3 .s and i; = -1/3. 

The plume theory of McCaffrey predicts a temperature gradient from 0.025 m to 1.5 m below the 
ceiling of only 2 °C, which compares well with the measured temperatures which were coincident 
within the uncertainty band. The predicted temperature at 0.025 m below the ceiling is 53 °C and at 
1.5 m below the ceiling is 55 °C. The predicted temperature at 0.025 m below the ceiling represents 
an underprediction of the excess temperature by 21 percent when compared to the measured value 
averaged over a 32 - 100 second time interval. 

Figure 10 shows the measured gas temperatures along the ceiling, 0.15 m (0.5 ft) below the ceiling at 
radial distances of 0, 1.5, 3.0, and 4.6 m (0, 5, 10, and 15 ft). It can be seen that there is no 
measured temperature gradient within the plume. 

4.2.2 Outside the fire test curtained area 

The measurement location at a radial distance of 12.5 m (41 ft) from the centerline of the fire, 
corresponding to a nondimensional radius of r/h = 0.41, is outside the plume region and on the other 
side of the fire test curtained area. None of the computer programs account for the flow around the 
draft curtains, thus this location can not be modeled accurately with present computer fire models. 

To model the draft curtain, and the flow around the curtain, a computational fluid dynamics model is 
necessary. 

4.2.3 Sensitivity to Heat Release Rate 

As was discussed in section 2.2 of this report, an analysis was conducted to determine the sensitivity 
of model predictions to the heat release input. Temperature predictions of the models were 
determined using a 20 percent higher heat release rate of 9900 kW (9380 Btu/s). Figure 1 1 shows 
the DETACT-QS, FPETOOL, and LAVENT predictions for the higher heat release rate at the 6. 1 m 
radial distance compared to the experimental data. This figure can be compared to figure 8 which 
shows the temperatures predicted for the same location with a heat release rate of 8250 kW. 
FPETOOL and DETACT-QS predicted a 9 percent (2 °C) higher excess temperature for the 20 
percent higher heat release rate. LAVENT predicted an 11 - 13 percent (2-4 ° C) increase in excess 
temperature. The higher heat release rate did change the temperature predictions; however, it did not 
change the overall results. That is, DETACT-QS and FPETOOL still underpredicted the gas 
temperatures. LAVENT predicted the temperatures well (closer to the experiment) at the position 
below the ceiling corresponding to the maximum temperature and velocity. 

4.3 Hot gas transport time 

The measured gas temperatures at the pan centerline, 1.5 m, 3.0 m, 4.6 m (5 ft, 10 ft, 15 ft) and out 
to 12.5 m (41 ft) to the east all show a significant jump in temperature at the 20 s data point (table 9 
givestime to first response of the thermocouples) that had not registered at the 16 s data point. This 
tells us that the smoke plume reached these points 18 ± 2 seconds after ignition. It is impossible to 
determine more exactly than the resolution of the data recorder. Also, the thermocouple 3.0 m (10 ft) 
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to the west of the pan centerline registered an increase in temperature 20 seconds after ignition. The 
thermocouple measurements combined with visual test observations, suggest the plume width to be 
approximately 16 m (52 ft). The calculated plume width (using Heskestad’s plume theory) at the 
ceiling is 15 m (49 ft) which agrees with this result. [19] 

Figure 12 shows the predicted gas temperature by LAVENT at a radial distance of 4.6 m (15 ft) and 
0. 15 m (0.5 ft) below the ceiling. This location was chosen because die measured temperature 
matched the predicted temperature so closely. Also shown in figure 12 is the same gas temperature 
prediction shifted 18 seconds to the right. The 18 second shift was selected to match the times of the 
predicted gas temperature rise to the times of the measured temperature rises. None of the computer 
fire models account for the transport time in the plume so changes in the input fire immediately affect 
the predicted temperatures at the ceiling. A similar shift could be applied to the FPETOOL and 
DETACT-QS results. 


5.0 COMPARISON OF FIELD MODEL (CFD) CALCULATIONS TO MEASUREMENTS 

The 30.4 m (100 ft) hangar curtained area was modelled using a heat source to simulate the pan fire. 
Computer speed effectively limits the total number of grid volumes available to model the problem. 

It was assumed that all the heat was released in the sixteen model grid volumes situated directly over 
the pan. Due to the size of the enclosure, only a portion of the curtained area was modelled. To 
reduce the number of grid volumes by a further factor of two, it was assumed that the fire was 
symmetric about a plane passing through the fire center and parallel to the draft curtains. 

Initial modelling efforts concentrated on the gas flow interacting with the draft curtains and assumed 
that the open trusses would have minimal impact on the gas flow. Heat flow to the ceiling was 
initially ignored due to the low expected temperatures at the ceiling. Later calculations included a 
conducting ceiling. Thermocouple measurements prior to the fire indicated that there was only a small 
thermal gradient present in the hangar. 

Several grid volume sizes were used in an effort to determine what effects, if any, were caused by the 
grid volume coarseness. A section view of the three dimensional grid used for some of the runs 
discussed later in the paper is displayed in figure 13. The grid structure is nonuniform in an attempt 
to m ax imize the grid volumes in regions where the temperature gradients are largest or where small 
structures are modelled. Details of the grid dimensions are given in table 7. 

Since the curtained area was not modelled in its entirety, pressure boundary conditions were used to 
bound the modelling regime at gas interfaces. At a pressure boundary condition, it is assumed that 
the inward or outward flow is continuous with upstream values of temperature, scalars, and 
turbulence used if flow is outward and ambient values used if flow is inward. The pressure 
boundaries were kept a minimum of twelve grid volumes away from the edge of the fire in order to 
minimize any error introduced by these boundaries. 

One of the major decisions in modelling the fire induced plume was the choice of turbulence model 
used in the calculation. For this work, the k-e model was used. This model has a number of user 
specified parameters which affect the ultimate viscosity and heat transfer in the plume. A number of 
investigators have proposed sets of parameters for the k-e model [20,21]. The default values used in 
the CFD model, HARWELL FLOW3D, for C, , Cj and C M are the parameters first proposed by 
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Launder and Spalding [17]. For buoyant plumes, an additional term, C^, to account for the 
production/destruction of turbulence due to buoyancy may be included. Additional changes have been 
proposed for some of the constants in order to match existing plume correlations with CFD 
modelling. Several sets of parameters for the k-c model, see table 6, were used in modelling the 
hangar fire. 

In modelling the fire for these calculations, it was assumed that 35 percent of the total heat release 
rate was radiated from the plume at the fire source. The r emaining 65 percent was then used as the 
convective heat release rate which provided the energy to drive the fire plume flow. The 35 percent 
radiative loss rate has been the standard value found in many computer fire models. For this fire, as 
mentioned earlier, the radiative loss rate may range from 20 to 35 percent [5]. 

The plume centerline temperature was determined in the following way. Unlike the thermocouple 
measurements fixed near the ceiling which maintained roughly constant temperatures throughout the 
first 150 seconds of the experiment, the thermocouples suspended at the geometrical plume center did 
not reach constant temperatures until about 100 seconds into the experiment. Visual inspection of the 
videotape suggests that the swaying of the thermocouple tree and some plume asymmetry may be 
responsible for the lack of constant temperatures early in the measurement. Therefore, the plume 
centerline temperature was obtained by averaging the thermocouple measurements taken between 100 
and 120 seconds which should give a "steady state" plume centerline temperature profile for the fully- 
involved pan fire. 

Figure 14 displays a comparison of the calculated and measured plume centerline temperatures as a 
function of height above die floor for each set of turbulence constants. The calculations are based on 
a time of 45 seconds after the start of the fire with time zero equal to 8250 kW (7810 Btu/s). The 
standard k-e constants with no buoyancy, k34, consistendy underpredict the centerline temperature 
and while the shape of the temperature profile in the lower plume is simulated, the upper plume is not 
well represented. Including the buoyancy term, k35, causes the model to overpredict the centerline 
temperatures and produces a much poorer agreement with height. Modifying two of the standard 
constants, calculations k36 and k38 provide better agreement with the temperature near the ceiling, 
but poor agreement with the overall shape of the centerline temperature. 

A comparison of the experimental data to the McCaffrey and Heskestad correlations [18,19] for 
plume centerline temperature are given in figure 15. The correlations, which are based on 
experimental measurements conducted in smaller spaces and axisymmetric conditions do not 
adequately represent the present data. The shape of the centerline profile does not provide for the flat 
temperature profile near the ceiling observed in the experiment. Calculations k34 - k38 are also 
compared with the McCaffrey result in figure 16. Calculations based on the turbulence constants 
used for k36 and k38 provide the best agreement with the McCaffrey correlation. 

One of the major difficulties in modelling the hangar is that the roof area is not only partitioned by 
draft curtains but also contains a number of beams, as mentioned earlier, which could act as 
additional smoke blockages and could account for at least part of the flat temperature profile observed 
near the ceiling. Figure 17 displays a comparison of a CFD calculation using the turbulence constants 
of k36 with the experimental data. For this calculation, the grid volumes were made smaller near the 
ceiling such that the 3.7 m (12 ft) draft curtain was described by 12 grid volumes rather than three as 
was done in the earlier calculations. The adiabatic assumption was relaxed and the ceiling was 
assumed to be conducting and made of steel which was assumed to be of thickness of 0.15 m (0.5 ft). 
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While the inside of the ceiling was corrugated steel, the makeup of the ceiling above the corrugation 
was not obtainable. Over the duration of the calculation, 30 s, the material above the steel was 
unimportant. It should be noted that the earlier calculations were run for 45 s while this calculation 
was run for 30 s. The shorter run time was used for the larger number of grid volumes after 
verifying, using the smaller number of grid volumes, that the calculation had essentially reached 
steady state at 30 s. Also included in the calculation were large I-beams which ran perpendicular to 
the draft curtain at the ceiling and effectively converted each ceiling compartment into a square. The 
additional framing trusses near the ceiling and the additional I-beams running perpendicular and 
parallel to the draft curtain at an elevation about 24 m (79 ft) above die floor were not modelled. 
These modifications provided an improvement to the calculation with the predicted temperature closer 
to the measured temperature near the ceiling and the deviation in shape away from the measured 
temperature being smaller at the lower elevations. 

The observation that both the plume correlations and the CFD computation suggest that the plume 
temperature should decrease with increasing height from the floor rather than remain a constant as 
was measured suggests that the experiment or the building geometry needs to be re-examined. While 
the draft curtains and ceiling beams may play some role in the temperature remaining constant near 
the ceiling, they should not provide much of an effect at heights equal to one-half the ceiling height 
for a building of this size at times early in the experiment. The explanation for the temperature 
variation with height may be experimental in nature. One possibility is that the plume is not 
maintainin g a completely vertical rise. In this case, thermocouples located at distances between 15 
and 25 m (49 and 82 ft) above the floor could be sampling parts of the plume that are cooler than at 
the centerline of the plume. Depending on how much the plume has swayed, a constant temperature 
profile with height could be experimentally produced using a single string of thermocouples located at 
the geometrical center of the fire as was done in this experiment. Based on plume correlations, the 
temperature half-width should be 1.9 m (6.2 ft) at the 15 m (49 ft) height above the floor. Over a 15 
m (49 ft) length, for the plume to sway 1 .9 m (6.2 ft) requires an angle of about 7 degrees as 
measured from the floor. Closer to the ceiling, the temperature half-width for the plume would be 
larger and the averaging effect would be less noticeable. Unfortunately, additional thermocouples 
were not available to check the plume symmetry at lower elevations. Visual inspection of the plume 
using the video playback suggests that the plume is fairly vertical but visual observations are line of 
sight dependent and particularly suspect when small angles are involved. Inspection of the flame in 
figure 5 suggests that there is some kinking of the flame which may continue upward until the plume 
impinges at the ceding. If the thermocouples are sampling across die plume diameter rather than the 
plume centerline, lower average temperatures will be obtained. 

A second possibility would be the sway of the thermocouple tree in response to the turbulence of the 
pl ume thus again sampling a lower average temperature than the maximum centerline temperature. 
The raw data shows a large variation of measured temperature with time indicating the unsteady 
nature of the measurement. Inspection of the videotapes indicate that the thermocouple tree was 
swaying in response to the fire turbulence. Whether these temperature variations were in part a 
response to the thermocouple tree swaying in the plume cannot be determined using the limited 
instrumentation associated with the plume. 

A third possibility would be a pre-existing temperature stratification which would slow the plume, 
inhibit entr ainm ent and possibly produce a constant temperature profile. The measured temperature 
stratification is at most 2 °C; hence, the plume should not have lost enough buoyancy to produce the 
observed experimental effect. 
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To summarize the results of this comparison, the CFD calculation using either k36 or k38 turbulence 
constants provide a fair representation of the plume centerline temperature at the ceiling. The more 
detailed grid volume CFD calculation provides a better representation of the plume centerline 
temperature near the ceiling than either the earlier CFD calculations or the two plume correlations. 

Temperature measurements along the ceiling between the plume center and the draft curtain were 
fairly constant in the experiment with the average temperature in the interval between 30 - 50 seconds 
only decreasing by about 4 °C for measurements 0.15 m (0.5 ft) beneath the ceiling. The CFD 
calculations at 39 seconds using the more detailed grid volumes at 0. 15 m (0.5 ft) beneath the ceiling 
gave a dip of 14 °C as shown in figure 18. At 12.6 m (41.3 ft), which is beyond the draft curtain, 
the temperature drops substantially for both the experimental measurement and the calculation with 
the calculated temperature in good agreement with the experiment at this location. 


6.0 RESULTS AND CONCLUSIONS 

1. Inside the plume region, within the curtained area, both the FPETOOL and DETACT-QS 
computer programs significantly underpredict the gas temperatures. The correlations of Heskestad 
and McCaffrey, and the calculation using the maximum temperature predicted by LAVENT also 
underpredicted the gas temperatures but were closer to the measured values. 

2. The CFD calculation using the modified k-e parameter provided the best agreement with 
experiment. 

3. For large spaces, a hot gas transport time delay should be incorporated into DETACT-QS, 
FPETOOL, and LAVENT. 

4. Established plume theory correlations such as those of McCaffrey for plume temperature and 
velocity in the r/h < 0.2 region need to be incorporated into the three computer fire models studied, 
since for high bay spaces the sprinkler or detector of interest would normally be located in this 
region. 

5. Reasonable agreement was found between the field model calculations using modified k-e 
constants and the measured centerline plume temperatures near to the ceiling. Further from the 
ceiling, the agreement between the calculations and the measurement was not as good. 

6. Reasonable agreement was found between the CFD calculation and temperature measurements 
along the ceiling. The temperature jump across the draft curtain can be qualitatively reproduced by 
the CFD calculation. 
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Table 1. 


Heat Release Rate as a Function of Time 

Time (s) 

Energy (kW) 

0 

0 

; 5 

4100 

8 

8250 

200 

8250 


Table 2. 


DETACT-QS Input Values 

Height of Ceiling Above Fuel (m) 

Distance of Detection from Axis of Fire (m) 

30.4 

0,1.5,3-0,4.6,6.1,12.5 

Initial Room Temperature (*C) 

22 

Detector Activation Temperature (* C) 

999 

Detector Response Time Index (RTI) (m-s) 1/2 

33.4 

Rate of Heat Release 

See table 1 
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Table 3. 


FPETOOL Input Values* 

ROOM 

1 SPACIAL DIMENSIONS 



Ceiling Height (m) 

30.4 


Room Length (m) 

81 


Room Width (m) 

12.5 

LINING MATERIALS (CEILING AND WALL) 



Material 

steel 


Thermal Conductance 



(kW/m-K) 

0.043 


Density (kg/m 2 ) 

7.85E+03 


Specific Heat 



(kJ/kg-K) 

0.48 


Thickness (mm) 

100 

SPRINKLER 


RTI (m/s) 1/2 

33.4 


Response Temperature (°C) 

999 ;..V 


Radial Distance (m) 

0,1.5,3.0,4.6,6.1,12.5 


SPRINKLER IS NOT A SIDEWALL 


FIRE 1 


Fire Height (m) 

0.15 


Heat of Combustion (kJ/kg) 

30,450 


Rate of Heat Release 

see table 1 

VENTS 



| DOOR/WINDOW 



Opening Height (m) 

26.8 


Opening Width (m) 



Sill Height 

0 


*No smoke detectors, HVAC, or program halts defined. 

1 Fire parameters listed are those differing from defaults program. 
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Table 4. 


LAVE NT Input Values 

ROOM PROPERTIES 

Ceiling Height ( m) 

30.4 

Room Length ( m) 


Room Width ( m) 

81 

Number of Vents 

0 

Curtain Length ( m) 

162 

Height for bottom 


of Curtains (m) 

26.8 

PHYSICAL PROPERTIES 

Ambient Temperature (K) 

295 

Material 

bare metal deck 

Material Thickness (m) 

0.1 

OUTPUT PARAMETERS 

Final Time (s) 

200 

Output Interval (s) 


FUSIBLE LINK PROPERTIES 

Radial Distance (m) 

0,1.5,30,4.6,6.1 

Distance Below Ceiling (m) 

0.15,0.30,0.91,1.5 

Radial Distance (m) 

12.5 

■ Distance Below Ceiling (m) 

0.025,1.2,2.4,3.7 

Radial Distance (m) 

0 

Distance Below Ceiling (m) 

0.025 

RTI (m-s) 1/2 

33.4 

Fuse Temperature (K) 

999 

|FIRE PROPERTIES 

Fire Height (m) 

0.15 

Fire Diameter (m) 


Fire Output as Function of Time 

see table 1 

SOLVER PROPERTIES 

program defaults 
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Table 5. 


Table 6. 

FLOW 3D INPUTS 

MODEL PHYSICS 

Turbulent Flow (k-e Model) 

Compressible Flow 
Heat Transfer 
Buoyant Flow 
Transient Flow 


GRID VOLUME DIMENSION 
21 x 28 x 90 grid volumes 


x: 

13 * 0.6858 m 

8 * 1.3716m 


■ 

y: 

10 * 1.3716 m 

8 * 0.6858 m 

10* 1.3716 

1 

z: 

40 * 0.1524 m 

30 * 0.6096 m 

20 * 0.3048 




AMBIENT VARIABLES 


Velocity (m/s) 
Gauge 


0 

Pressure (Pa) 


0 

Density (kg/m 3 ) 


1.2 

Viscosity (kg/m-s) 


1.91E-05 

K (J) 


1.00E-10 

« a/s) 


1.00E-11 

Temperature (K) 
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BOUNDARY CONDITIONS 

Pressure Boundary (Sides) 


Solid (floor and non-conducting ceilings) 


Conducting Solid (for conducting ceilings) 


































Table 7. 


Plume Symmetry 


Mean Temp °C (32-200 s) 

Standard Deviation 

Location 

Depth Below 
Ceiling (m) 

Test 1 

Test 2 

Test 3 

Test 1 

Test 2 

Test 3 

East 3.0 m 

0.15 

63 

59 

64 

7 

6 

7 

from pan 




61 

Mlimi 


7 

centerline 

0.91 

63 

61 

65 

6 

4 

6 




wmmm 

6 5; 



7 

West 3.0 m 

0.15 

53 

59 

57 

4 

6 

B 

from pan 

0.3 


mimm 

^ . 54 ,/: 

: ■■■ ■ -;';;4' H 

:^6;= . - 

B 

centerline 

0.91 

52 

58 

57 

4 

5 

B 
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Table 8. 


Mean Gas Temperatures (32-100 s) 


Depth Below 

T 

x measured 

Standard 

Location 

Ceiling (m) 

(*C) 

Deviation 

Pan 




Centerline 




0 m 


64 

(6) 

East 1 .5 m 

0.15 

61 

(6) 

(r/h = 0.05) 


59 

(6) 


0.91 

63 

(5) 



62 

( 5 ) 

East 3.0 m 

0.15 

59 

(6) 

(r/h = 0.1) 


55 

m 


0.91 

61 

(4) 


1.5 

61 

( 5 ) 

East 4.6 m 

0.15 

57 


(r/h = 0.15) 

0.3 

55 



0.91 

59 

El 1 


1.5 

57 

E9 1 

East 6. 1 m 

0.15 

59 

(4) 

(r/h = 0.20) 

0.3 

59v: 

(4) 


0.91 

64 

(6) 


,'.1-5 

63 

(5) 

East 12.5 m 

0.025 

33 

(4) 

(r/h = 0.41) 

1.2 

35 

(3) 


2.4 

42 

(7) 


3.7 


(5) 
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Table 9. 


Hot Gas Transport Time 

Radial Distance 
From Pan Centerline (m) 

Time to Thermocouple 
Response (s) 

0 

20 


20 

East 3.0 

20 

East 4.6 

20 

East 6.1 

20 

East 6.2 (under draft Curtain) 

20 

East 12.5 

East 18.9 (under draft curtain) 

20 



West 3.0 

20 

West 5.5 

24 
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igure 2. Plan view of hangar 
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Figure 4. Fire test pan 




33 


Figure 5. Picture of fire at steady state 
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Figure 6. Instrumentation layout 



Temperature (°C) 



Time (s) 


Figure 7. Measured vs. predicted gas temperatures at a radial distance of 1.5 m from 
the fire centerline, r/h = 0.05. d = distance below the ceiling 
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Figure 8. Measured vs. predicted gas temperatures at a radial distance of 6.1 m from 
the fire centerline, r/h = 0.20. d = distance below the ceiling 
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Figure 9. Illustration of plume and ceiling jet regions of a typical fire 
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Temperature (°C) 



Time (s) 


Figure 10. Temperature distribution within the plume measured at radial distances (r) of 
0, 1.5, 3.0, 4.6, and 6.1 m from the fire centerline, 0.15 m below the 
ceiling.. 


38 




Temperature (°C) 
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40 
20 
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Figure 11. Measured vs. predicted gas temperatures at a radial distance of 6.1 m from 
fire centerline. Temperatures predicted using a 20 percent higher heat release 
rate of 9900 kW 
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Temperature (°C) 



Temperature at 4.6 m radial distance from fire centerline with time shift, 
r/h = 0.15 d = distance below the ceiling 



z 

A 


Draft curtain 


Fire 



19.9 m 


T 


* 


Figure 13. Two dimensional section of a three dimensional grid used for the CFD calculations 
representing a portion of the fire test curtained area and the east adjacent curtained 
area 
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30.5 m 



Figure 14. Comparison of CFD calculations using four different sets of k-c constants to 
the measured plume centerline temperatures 
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McCaffrey plume correlations to the measured 



Temperature (°C) 



10 15 20 25 30 

Height Above Floor (m) 


Figure 16 . Comparison of calculations using four different sets of k-e constants to the 
plume correlation of McCaffrey 



T emperature (°C) 
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Figure 17. 


Comparison of a CFD calculation using the turbulence constants of k36 and a 
more detailed grid with the measured plume centerline temperatures 


I — X — i 






Figure 18. Comparison of the CFD calculated ceiling temperatures using a detailed grid 
and the k36 turbulence constants to the measured temperatures at various 
radial distances from the fire centerline and 0. IS m below the ceiling. 
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1.0 INTRODUCTION 


NASA has numerous large, high bay spaces that are used to perform a variety of different functions. 
Many of these functions are critical to NASA operations. Examples of high bay spaces at NASA 
include those used for clean rooms, shuttle simulators, assembly/storage, vacuum and vibration 
chambers, vehicle assembly, and/or testing facilities with payloads. These spaces represent some of 
the most difficult fire protection challenges in that detection of a fire in a large space may be delayed 
due to the distance smoke and products of combustion must travel to reach the detector, the large 
amount of ambient air for smoke dilution, the high dollar value of these spaces and the low damage 
threshold of a clean room. Some of these spaces also involve forced airflow. 

As was shown in section A of this report, current computer fire models do not provide accurate 
temperature predictions in a high bay space where the region of interest is within the plume. Design 
for fire protection in these spaces at present relies upon "best estimate" guess work. Therefore, it is 
important for NASA to develop a methodology by which to determine fire protection options for these 
spaces. The methodology must be in a suitable form to be incorporated into the NASA fire protection 
manual and be easily usable by the field engineers. This project represents the first step toward 
developing such a methodology. An existing NASA high bay space was modelled using a field model 
to determine how the existing fire detection devices would respond to various size fires in the space. 

The CFD model, HARWELL FLOW3D [1], used to model the NASA high bay space has already 
been successfully applied to the analysis of fire problems in both forced and unforced flow situations 
at NIST. Compartment fires with single and three room configurations were successfully modeled [2] 
as part of a validation project, section A of this report extended the validation efforts to high bay 
spaces. The model was used to successfully simulate an experiment [3] used to determine the impact 
on the ceiling jet flow of beams on a flat ceiling. Finally, die model has been used to analyze and 
solve a fire problem associated with the Navy’s 19F4 Fire Fighting Trainer where winds blowing in 
from the ocean adversely impacted the operation of the trainer [4], 

A clean room at Goddard Space Flight Center in Greenbelt, MD, was chosen to be analyzed for 
detector placement based on the premise that the forced air flow through the room presented the 
possibility of adversely effecting the detection of smoke from a fire. Because of the air flow in the 
room, models such as DETACT-QS, FPETOOL, and LAVENT [5-7] used for the earlier analysis are 
not appropriate for use in this space. These computer fire models assume that the plumes are 
axisymmetric, neglect flow velocities in both the upper and lower layers and make the assumption that 
an upper layer is formed regardless of fire size and flow conditions. The air flow in this clean room 
may have sufficient velocity to negate all of these assumptions for small fires. The purpose of this 
analysis will be to study the impact of fire size on the distribution of hot gas and detector activation in 
the clean room. 

2.0 NASA GODDARD HIGH BAY CLEAN ROOM 

The Goddard clean room is the largest known laminar flow clean room in the world. The space 
measures 38.1 x 30.5 x 27.4 meters (125 x 100 x 90 ft) in dimension. The room is rated as a class 
10 clean room and is large enough to accommodate two full size shuttle payloads simultaneously. 

Two decks are located on the side of the room opposite the air intake wall. Two 35 ton cranes with 
hook heights of 21 and 24.4 m (69 and 80 ft) are mounted on beams near the ceiling. HEPA filters 
are used to remove particulate matter at the air intake wall. A perforated return air grille is used for 


49 



the outlet wall which is located opposite the intake wall. The laminar air flow in the room is uniform 
from floor to ceiling in a horizontal direction and is maintained at a rate of 0.508 m/s (1.67 ft/s) by a 
computer operated mechanical control system. The concrete floor is covered by a polyurethene 
coating. The clean room is used for assembly and testing of space flight hardware. 

A diagram of the room showing the gridding is given in figure 1. The solid decks are shown as 
rectangles extending out from the right side of the grid. The cranes are approximated by the two 
squares inserted in the upper right portion of the grid. Figure 1 represents one vertical plane 
extending from the floor to the ceiling of the clean room. The grid is actually a three dimensional 
grid of size 25 x 20 x 24. The large grid volumes are cubic with each side having a length of 1 .524 
m (5 ft) and the smaller grid volumes have a vertical length of 0.762 m (2.5 ft) with a square base of 
1 .524 m (5 ft) on a side. 

The smoke detectors on the outlet wall are spaced every 3.66 m (12 ft) in both the horizontal and 
vertical direction. Since the two decks essentially compartmentalize the outlet wall into three spaces, 
the initial detector design required the placement of at least one row of detectors in each space. 

There are no detectors on the ceiling. 

3.0 ANALYSIS 

3.1 FIRE WITH FORCED AIRFLOW 

Smoke flow from a fire will be affected by both the fire size, which will impact buoyancy, and the 
forced airflow in the room. Large fires, which will create strong buoyant forces, should produce a 
plume which will propel most of the smoke toward the ceiling and form a ceiling layer. Small fires 
will have weak buoyant forces which may result in very little upward movement before the ambient 
air flow causes the smoke to exit the room. The horizontal smoke spread perpendicular to the flow 
velocity, may be reduced with a small fire, hence detector spacing for small fires may be more 
critical than for large fires. 

In the model, the fire was placed three grids off of the floor and in the far center of the room near 
the inlet wall. The location of the fire off the floor was suggested by the raised assembly areas which 
were in active use when the clean room was visited. The fire was modelled by uniformly distributing 
the heat release rate in four adjacent grid volumes. Due to the geometrical simplicity of the room, it 
was felt that varying the location of the fire was not necessary. The default values for the k-e 
turbulence model were used in these calculations. 

The model calculations give the location of the heated gas but do not explicitly track smoke particles. 
A connection between the calculated hot gas distribution and the smoke distribution is therefore 
required in order to use the hot gas distribution to track the smoke particles. It is assumed in this 
analysis that the smoke distribution is coincident with the distribution of hot gas. This assumption is 
supported by the experiments of Heskestad and Delichatsios [8] which showed that for a fire under a 
flat smooth ceiling, the optical density (smoke density) is correlated with the local rise in temperature. 
While there was no forced airflow present in these experiments, the smoke did follow the gas flow 
produced by the buoyant force of the fire. Unless there is substantial deposition of smoke it would be 
expected that in forced flow situations, the smoke would again follow the gas flow and that the smoke 
density would correlate with the local temperature rise. Therefore, for the analysis used here, it was 
assumed that the excess temperature plots were sufficient to locate the smoke particles. 
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Three fire sizes, 32 MW, 400 kW, and 40 kW (30000 Btu/s, 380 Btu/s, 38 Btu/s) were used to 
investigate the impact of fire size on detectability. The 32 MW (30000 Btu/s) fire was chosen as it 
represents the maximum heat release rate for the default fire in the model, LAVENT. The other fire 
sizes represent smaller, more likely fires for the clean room with the 40 kW (38 Btu/s) fire 
approximating the size fire found in a wastepaper basket. Contour plots, in a vertical plane passing 
through the fire, of excess temperature (temperature above ambient) are given in figures 2, 3, and 4 
for the three fires. The darker shaded areas represent higher temperature areas which should contain 
higher smoke concentrations. Clear white areas will contain no smoke. The small values for the 
excess temperatures are due to the large grid volumes used to model these very small fires in this 
large space. 

From the figures, it can be seen that the larger two fires produce sufficient buoyancy to allow the 
heated gas to reach the ceiling and hence may be detected by at least some of the wall mounted 
detectors near the ceiling. For the largest of the fires, there is potential for plume driven circulation 
patterns to set up where some of the smoke detectors at lower elevations may not detect the fire. 

The small fire did not possess enough buoyancy for the fire plume to rise to the ceiling by the time 
the heated gas exited the room via the outlet vents. Hence, the detection of a small fire in this space 
will require wall mounted detectors which are approximately at the height of the fire. The spread of 
the heated gas in the horizontal direction at the base of the fire is shown in figure 5. There appears 
to be adequate horizontal spread for the smoke to encounter at least two smoke detectors. 

The computer fire models used in section A assume that the fire plume rises vertically to the ceiling. 
The forced ventilation effects on the fire plume predicted by the CFD calculations for small fires 
would not be predicted by these fire models. 

While the distribution of the smoke detectors may be adequate, the question still remains concerning 
whether the sensitivity of the detectors are adequate for small fires. The wall-mounted detectors in 
the clean room are photoelectric smoke detectors which have a rating of 5.2 - 8.1 percent obscuration 
per meter (1 .9 - 4.8 percent obscuration per foot). This rating is related to the optical density per 
meter of the smoke by (9] 


Where I x is the observed intensity, I oX is the incident intensity, D is the optical density per meter, and 
L is the path length. D is related to the fire by the relation 

D=D m (AmlVJ 


Where D m is the mass optical density. Am is the mass consumed in the fire, and V c is the smoke 
volume. 

The above two equations must be evaluated for a given fire in order to estimate the optical density at 
the detectors. The values of D m and Am depend on the fuel and also on the heat release rate. Fuels 
typically have heats of combustion between 16 and 40 kJ/g and produce mass optical densities 
between 0.12 and 1.4 m 2 /g (590 and 6800 ft 2 /lb). Fuels found in the clean room may consist of 
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plastics, mylars, metal, and wood. For the calculations given below, it will be assumed that the 
material heat of combustion will be 30 kJ/g (28000 Btu/lb) and the mass optical density produced will 
be 0.5 m 2 /g (2400 ft 2 /lb). This will give a middle range result with some materials producing more 
smoke and others less smoke. The smoke mass per volume Am/V c will be estimated using the 
following relation 


Am Q 

K H A V 


where H c is the heat of combustion, 0 is the heat release rate, A is the area occupied by the smoke 
perpendicular to the forced airflow, and v is the velocity of the forced airflow in the clean room. The 
area A at the detector wall is determined from the field model calculations by assuming that the 
smoke is essentially contained in the region at the wall where the excess temperature is greater than 
1/3 of the maximum excess temperature. 

The smoke from the 40 kW (38 Btu/s) fire impacts an area on the outlet wall of approximately 140 
m 2 (1500 ft 2 ), approximately 2/3 of the area under the lower deck. Using the above mass optical 
density and heat release requires a detector with a rating of 2.3 percent obscuration per meter (0.7 
percent per ft) which is significantly below the sensitivity of the existing detectors. The 400 kW (380 
Btu/s) and 32 MW (30000 Btu/s) fires had perpendicular areas which approximately equaled the clean 
room wall, 840 m 2 (9000 ft 2 ). The 400 kW (380 Btu/s) fire would require a rating of 3.9 percent per 
meter (1.2 percent per ft) or 75 percent of the existing detector rating while the 32 MW fire required 
a rating of 200 percent obscuration per meter (61 percent per ft). 

Both the 40 kW (38 Btu/s) and the 400 kW (380 Btu/s) fires generated estimated smoke 
concentrations below the detection limit of the existing smoke detectors in the clean room. This 
conclusion depends on the smoke following the heated gas and being well mixed. 

3.2 FIRE WITHOUT FORCED AIRFLOW 

The 400 kW (380 Btu/s) fire was used to demonstrate how the analysis can change when the forced 
airflow is absent in the clean room. The heat source representing the fire was positioned as in the 
earlier analysis but a smaller grid was used to define the plume more accurately. Two openings, 
simulating open doors, centered at the floor on opposite walls are used to vent the room. Figures 6 
and 7 present the calculated temperature distribution at times of 50 and 100 s after the fire ignition. 

As can be seen from the figures, a vertical plume rises from the fire source, is deflected by the 
ceiling and forms a ceiling jet along the ceiling. Since the fire is not symmetric in the room, the 
ceiling jet strikes the near wall first and sets up a circulation pattern there. The smoke is 
concentrated in the ceiling jet and does not diffuse as it does in the forced airflow environment. The 
gas temperature without forced airflow, see figures 3 and 6, is substantially higher than is calculated 
in the forced airflow simulation. In these situations, smoke detectors located at substantial distances 
beneath the ceiling would be of little use in the early detection of the fire. Because the smoke is 
concentrated in the ceiling jet, the requirements on smoke detector sensitivity near the ceiling would 
be less severe than in the forced airflow environment. 
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4.0 RESULTS AND CONCLUSIONS 


1 . Computer fire models examined in section A of this report are inappropriate to use for detector 
analysis in forced air environments. 

2. The design fire size is of extreme importance in structures where there is a forced air flow. For 
the NASA Goddard clean room, small fire plumes would not reach the ceiling, hence the distribution 
and sensitivity of side wall detectors becomes important. 

3. Large fires in structures containing forced air flow may provide sufficient buoyancy for the fire 
plume to reach the ceiling. In these situations, wall mounted detectors may be shielded by 
obstructions from the smoke and not detect the fire. 

4. Detector sensitivity in the NASA Goddard clean room should be increased for small fire detection. 
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Figure 1. Vertical slice of a typical grid used for the NASA clean room calculation 
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Figure 2. Temperature contours on a vertical plane through the fire for the 32 MW fire 
calculation 
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Figure 4. Temperature contours on a vertical plane through the fire for the 40 kW fire 
calculation 
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Figure 5. Temperature contours on a horizontal plane centered on the grid directly above the 40 
kW fire 
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Figure 6. 


Temperature contours on a 
400 kW fire 
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